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Abstract 

In this paper we investigate the scalar mode of first-order metric perturbations over spatially 
flat FRW spacetime when the holonomy correction is taken into account in the semi-classical 
framework of loop quantum cosmology. By means of the Hamiltonian derivation, the cosmological 
perturbation equations is obtained in longitudinal gauge. It turns out that in the presence of metric 
perturbation the holonomy effects influence both background and perturbations, and contribute 
the non-trivial terms Shi and Sh2 hi the cosmological perturbation equations. 
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I. INTRODUCTION 



In loop quantum cosmology two main quantum gravity effects lead to remarkable modifi- 
cations to the standard description of the early universe(for a detailed review, see Ref. 
One is due to the holonomy correction and the other is due to the inverse volume correction. 
Such modifications can successfully avoid the Big Bang singularity {2 5|, and replace it by 



the Big Bounce even at the semi-classical level p, Ql- In addition, it is very interesting to 
notice that quantum gravity effects may lead to the occurrence of the super-inflationary 
phase js]. As shown in Ref. jsl, such a super-inflationary phase can also resolve the horizon 
problem with only a few number of e-foldings. Therefore, it is possible to construct a phase 
of inflation or an alternative to inflation in the framework of loop quantum cosmology. 

As we all known, the inflationary phase is crucial for understanding the structure forma- 
tion and anisotropics of the CMB. In order to address these issues in the framework of loop 
quantum cosmology, we must consider the cosmological perturbation theory with modifica- 
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by means 



tions due to quantum gravity effects. In the earlier work by Bojowald et al. 
of the Hamiltonian derivation they have obtained the cosmological perturbation equation 
with inverse volume corrections for scalar modes in longitudinal gauge. They show that 
super-horizon curvature perturbations are not preserved. Recently, they have also derived 
the gauge-invariant quantities and the corresponding gauge-invariant cosmological pertur- 
bation equations with inverse volume corrections for scalar modes 
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13] . In addition. 



the vector modes and tensor modes with corrections from loop quantum gravity have been 

nn 

investigated [1^ 115[ . 

At the same time, some pioneer work have already been devoted to understanding the 
— p. powe. spect^un, in the perturbation theory of LQC fl First of ah, rn 

Ref. |9|, |18[, it is shown that a scale invariant spectrum can be obtained. More importantly 
these attempts imply that the quantum gravity effects may leave an imprints on the power 
spectrum which can be potentially detected in the future experiments such as the Planck 
satellite. However, above considerations are restricted to the scalar field perturbations with 
fixed background. To provide a complete and more precise understanding on the perturba- 
tion theory in loop cosmology, it is essential to take the metric perturbation into account. 
Along this direction it is worthwhile to point out that another potential observables, pri- 
mordial gravitational waves have already been investigated intensively in LQC 23 |. 
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Although, in Ref. |10l-ll3| the cosmological perturbation equations with inverse volume 
corrections have been derived in longitudinal gauge and gauge-invariant manner respectively, 
the metric perturbations with holonomy corrections is still absent. In the present paper, by 
means of the Hamiltonian derivation, we will derive the cosmological perturbation equations 
with holonomy corrections in longitudinal gauge. 

The outline of our paper is the following. For comparison, we firstly present a brief 
review on the perturbation equations in standard classical cosmology in section 11. After 
introducing the basic variables in loop cosmology in section III, we will demonstrate a 
detailed derivation on the cosmological perturbation equation with holonomy corrections in 
section IV. The discussion is given in section VI. 

II. THE CLASSICAL COSMOLOGICAL PERTURBATION EQUATIONS 

Before proceeding to the effective loop quantum cosmology with holonomy corrections, 
we first briefly review the classical perturbation equations in standard cosmology. A detailed 



24| . Let us now consider a spatially fiat background metric 



derivation can be found in Ref. 
of FRW type 

(is2 _ a\r]){-dr]'^ + dabdx^dx^) . (1) 

where r] is the conformal time. The spatial part of the metric describes isotropic and homo- 
geneous 3-surfaces. Then one can perturb the background metric 

ds"^ = a\r]) [-(1 + 2<!>)dr]^ + (1 - 2^)5«fedx"dx*] . (2) 

Here we only consider the scalar modes in longitudinal gauge, which is thus diagonal. 
Through this paper, we will consider the scalar field ip as the matter source. Expanding the 
Einstein's equation linearly, one can obtain the cosmological perturbation equation 

_ 3e$ - (H + 2e2)$ = AnGi^Sip + pV^{ip)5ip) , (3) 

$ + 3e<i) + (e + 2lf)(^ = 4TTG{<^6ip - pV^{p)6ip) , (4) 
dai^ + M^)=AnG^Sip,a , (5) 

where a dot denotes a derivative with respect to the conformal time r]. H is the Hubble 
expansion rate in the conformal time, and for later convenience, we have identified with 
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p which is introduced in f llOp . Note that in the case of vanishing anisotropic stresses, two 
scalar functions $ and coincide, $ = Therefore, in above equations we have set $ = 
which simphfies the equations considerably^. Moreover, among these equations above only 
two of them are independent. Combining these equations, one can obtain the following 
second order differential equation for $ 

$ _ + (6H + 2p^^^4^)$ + (2e + AE? + ^5:^4^H)<I> = . (6) 

LP Cp 

In addition, the background and the perturbed Klein-Gordon equation can respectively 
expressed as 

(^ + 2H(^ + p\/^(^) = , (7) 
Sip + 2m6ip - yHip + pV^^{^) + 2pV^{ip)^ - 4(^<i> = . (8) 

III. THE BASIC VARIABLES 

Now we intend to study the scalar mode of first-order metric perturbations around spa- 
tially fiat FRW spacetime when the holonomy corrections is taken into account in the semi- 
classical framework of loop quantum cosmology. To derive the cosmological perturbation 



equations we adopt the Hamiltonian approach which has 



Deen developed in the effective 
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12|. We summarize the basic 



loop quantum cosmology with inverse triad corrections 
idea and steps as follows. 

In loop quantum gravity, instead of the spatial metric qab-, a densitized triad is pri- 
marily used, which satisfies E^E\ = q"-^detq. Moreover, in the canonical formulation the 
space-time metric is given by 

ds"^ = -N^drf + qab^dx" + N''dr]){dx'' + N^dT]) , (9) 

where and A^" are lapse function and shift vector respectively. 

By comparing the above equation with the FRW metric ([T]), the background variables, 
N, N°- and E^, can be expressed as respectively 

N = ^-N- = Q-E'^=p5t . (10) 



^ However we must point out that, as a matter of fact, $ = ^' is a consequence of equations of motion, 
which can also be seen in this paper 
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where the background variables are denoted with a bar, which describe smoothed out, spatial 
averaged quantities. Another background variable, the extrinsic curvature components K^, 
can be derived from the following relation 

Kab = -^{qab - 2D^,N^)) = aSab . (11) 

where D is the covariant spatial derivation. Thus, the extrinsic curvature can be expressed 

as 

K = -^=Kab = =■■ H . (12) 



\det{E: 



2p " 



In equation (1121) . we have defined the background extrinsic curvature as k =: 



which can also be obtained from the background equations of motion [12|. Therefore, in 
classical FRW background, the extrinsic curvature is nothing but the conformal Hubble 
parameter H. However, in the effective loop quantum cosmology, the relation between the 
extrinsic curvature and the conformal Hubble parameter will change due to quantum gravity 
corrections, which we will see in the next section. 

The canonical perturbed variables can be related to the perturbed metric variables by 
comparing the perturbed metric ([2]) with the canonical one([9]). It turns out that the per- 
turbed triad is given by 

6E^ = -2p^6f , (13) 
and the perturbed lapse function is 

5N = iV$. (14) 

As shown in the above, the extrinsic curvature components can be diagonal, thus it can 
be expanded as 

K = K: + 5K = k5: + 5K:. (15) 

The perturbed extrinsic curvature will be derived from the equation of motion in the 
following. We can assume that SEf and do not have homogeneous modes, namely 

[ SE^Sid^x = 0, [ Sirjfd^x = . (16) 



And the homogeneous mode is defined by 



I EtS'j'x, k=^ I KlSfd'x , (17) 



E 



where we integrate over a bounded region of coordinate size Vo = J.^ d?x. Then we can 
construct the Poisson brackets of the background and perturbed variables 

SttG 



13|, 



(18) 



In addition, we point out that the similar conditions will be required in the perturbed 
lapse 6N, the scalar field 6lp and conjugate momentum Sir such that 

6N£x = 0, / 6<^d^x = 0, / 67t£x = , (19) 

which is used in expanding the Hamiltonian constraint. While the homogeneous mode of 
the scalar field and its conjugate momentum is 



ip = — f Lpcfx, vr = — f Trd^x . 
H Jt. K) is 



(20) 



Therefore, the Poisson brackets of the background and perturbed variables of scalar field 



IS 



vr} = — , {S^{x), S7T{y)} = 5^{x~y) . 

6Vo 



(21) 



IV. THE COSMOLOGICAL PERTURBATION THEORY WITH HOLONOMY 
CORRECTIONS 



Now we turn to the derivation of the cosmological perturbation theory in the effective 
loop quantum cosmology with holonomy corrections. For more details on the Hamilton 
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ij. 



cosmological perturbation theory, we refer to Ref. 

Thanks to the holonomy corrections, in the isotropic and homogeneous models, the ef- 
fective Hamiltonian can be obtained at the phenomenological level by simply replacing the 
background Ashtekar connection k by where 7 is the Barbero-Immirzi parameter. 



may be a function of p. More 



25|. 



The parameter Ji depends on the quantization scheme anc 
discussions about the parameter fi, we can refer to Ref. 

However, when the inhomogeneities are taken into account, it is no longer true. To study 
the effects of holonomy corrections on inhomogeneous perturbations, we similarly substitute 
the appearance of k in the classical Hamiltonian by a general form where m is an 

I 1 I 1 

integer. In the context of vector modes [ij] and tensor modes 15|, due to the requirement 



of the anomaly cancellation, we can fix the parameter m. Since the evolution of all modes 



should be generated by one general Hamiltonian constraint, it would be reasonable to use the 
values found for vector modes and tensor modes also for scalars. However, it must been also 
pointed out that the restrictions of anomaly cancellation from the vector modes and tensor 
modes has not been checked for scalars in the presence of holonomy corrections. Complete 
consistency is realized only if all modes can be anomaly-free with holonomy corrections for 
specific parameter values. 

Subsequently, one can write down the expressions for the gravitational Hamiltonian den- 
sity in a similar manner "H^ = 7/^°^ + T-Lq^^ + "H^^^ with 

1 ,sin^7^^,2,_,,,,2 ^'So xz7c^,o .j;.^^ 



+ ^^^:^n^my - ^.(d^m^Et) , (22) 

where the superscript "/i" represents the holonomy corrections and the corresponding clas- 
sical expressions can be found in Ref. 12, isl. We now only consider the scalar field as the 



matter source. Its Hamiltonian density expands as T-Lm = 'Hfj + "H^^ + T-L^m ■ Since the 
matter is free from the holonomy corrections, the expressions of scalar field Hamiltonian 
density, Hm = 'Hn + "Hv + "H^, expanding up to the second order, are as the classical cases 
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13|, 



and 



^1°) = ^2^' = 0, = P'^'vm, (23) 



-1/(2) ^ l^^TT^ n6nSi6E^ 1 ( {SjdE^)' 5%5E'^^5E f 

= \^/p5'''^a5^^,5v , 
= ^P^/V,,(^)V+P^/V,(^)5^^ 



-P"'V{^) I - - ' I . (25) 
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A. The background equations 



In the isotropic and homogeneous FRW background, the diffeomorphism constraint van- 
ishes. Therefore background equations are generated only by the background Hamiltonian 
constraint, which can be expressed as 

H'^i^^N] = -i- / d^xNinP + i%TiG{nf + nf)] . (26) 

iOTTG Jj^ 

Thus the exphcit expression for the background Hamiltonian constraint is 

3 sh^ *^ (27) 

SvrG /i7 2p-^i^ 

Then, by means of Poisson bracket, we can derive the equation of motion for the gravi- 
tational variables k and p. 

I = {K H^lN]} = _[1 + ^,G[-f^+pVm . (28) 

2 /i7 op fi^ 2p^ 

p = {p,H'^'m} = 2pC^^^). (29) 

Similarly, the equation of motion for scalar field (f and its conjugate momentum field n 
can also be derived as 

^ = {^,H'^(^\N]} = - . (30) 

P 

^ = {7r,//'^W[iV]} = -p%(^) . (31) 

Note that in above Poisson brackets, we have used the relation N = y/p. Substituting 
the relation fl30|l into the constraint equation fl27|) gives rise to the corrected Friedmann 
equation 

:V +pyw] ■ (32) 



pl 3 2 

At the same time, equation (l28l) is just the corrected Raychaudhuri equation 

*4(^)^+4(^)^-4.G|4.m^)l . (33) 

In the classical limit, /i — > 0, above two equations can be reduced to the Friedmann and 
Raychaudhuri equation in the standard cosmology. Finally, the Klein-Gordon equation can 
be derived from Eqs. (|30|). (|3T1) and (|29|) 

^ + 2(— ^)^ + pV^{^) = . (34) 
2/17 
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In addition, from the equation of motion (129|) . one can find that the extrinsic curvature 
k is related to the conformal Hubble parameter EI by 

'J^^ = l=:m. (35) 

Therefore, due to the holonomy corrections, the conformal Hubble parameter EI is not 
simply equal to the extrinsic curvature k but receives corrections. For consistency, in our 
next derivation we will continuously use the extrinsic curvature k rather than the conformal 
Hubble parameter. Only at the end, we will use the conformal Hubble parameter EI instead 
of '^"i-^^^ in the perturbation equations. 



B. The perturbed equations 

In this subsection, we will derive the cosmological perturbation equation with holonomy 
corrections. Firstly we will derive the equations of motion of perturbed variables. In the 
canonical formulation, the equation of motion of any phase space function / is determined 
by Poisson bracket, / = {f,H}. Here H is the total Hamiltonian, which is a sum of the 
Hamiltonian constraint H[N] and the diffeomorphism constraint D[N"-], H = H[N]+D[N"-]. 
Since the zero-order and first-order shift vectors vanish, the diffeomorphism constraints is 
identically satisfied up to the second-order. Thus, the equations of motion of the perturbed 
variables are only generated by the Hamiltonian constraint. The perturbed Hamiltonian 
constraint up to the second-order is written as H'^[N] = H'^[N] + H'^ISN] with 

^'[^^] = TTTT^ / d'xdNin^^^ + 167rG(?/« . (36) 

iOTTCjr Jy, 

Note that we have used the conditions that the perturbed variables do not have homo- 
geneous modes as described in Eq. flT6|) and f|T9l) . As well, we input the boundary condition 
requiring that the integration over the boundary vanishes, namely 

/ iV[H^i + IQnG^Hl + K)] = 0, / 6N[H''^ + lQixG{Hl + K)^ = ■ (37) 

Therefore, the equations of motion of perturbed variables are generated only by the 
second order part of Hamiltonian constraints. Thus, we can arrive at the equation of motion 
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of the perturbed variables by means of the Poisson bracket 



5K = {SKIH^[N]+H''[6N]} 

N . , sin 2/X7/C . ^ . 1 ,sin/i7A;,2„^rft.i,„- 1 /Sin/i7/c,2/,^rf„i.- „• ^^j, 

2 \/p VP 

+pV{v){\6Et6',6l - 5Et6%)] + AixG5N[-\^^ + Vpvm^l (38) 



5E^ = {6E^,H''[N] + H^ldN]} 

= ^[-p5Ki5^5] + p{SKiS^)6t + C^^^)^En + 25NC^^^)VpSt, (39) 

^ {5^, H'[N] + H'[SN]} = ^,{Sn- *^) + , (40) 

Sn = {Sn, H^[N] + H^[5N]} = ^[p^vH^ - p^V^^6^ - ^p^V^^5E;;5i] . (41) 



Furthermore, using Eqs. (jT3l) and ( 1291) . we can obtain the perturbed extrinsic curvature 
6K^^ from equation ( l39i) . 

sK = ~6:[^ + C^^^){^ + n- (42) 

2/i7 

Similarly, using Eq. (fT3|l . equations ( 140|) and ( 14T|) can be respectively reexpressed as 

^ + ![(3^ + $) . (43) 
p p 

Sn = pVH^ - p'^V^^Sif + 3pV^^ . (44) 

Now, we derive the Hamiltonian's equation using the equation of motion of SK^^. Col- 
lecting the expressions SE^ ([13]), 6Ki (02]) , 5N(^, and equations ([30]), ([H, one can 
obtain 

+ ^_ ^ 2^ + $ + cos2/i7A; - - A: + ' f + -{kcos2fL-fk 

2/^7 2 2/^7 fj, 2/i7 

_i _ a sm^ + $ ) +pvm^ - ^)}s: + a,a^(<i. - ^) 

2 ajo /i7 

= AivGi^dip - pV^i^)Sip) . (45) 
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When deriving this equation, we have used the relation 



AnGip = (— ) -k-p—{— ) , (46) 

fi^ op /i7 



which can be obtained from the corrected Friedmann equation fl32|) and Raychaudhuri equa- 
tion (133|) . From equation (145|) . we can read the off-diagonal equation 

a„a*[$-^]=0, (47) 

which implies $ = Therefore, in the following derivation, we will identify $ with 
Then the diagonal equation gives 

, sin 2/17 A;,- r/ r -.t ,sin2/i7A;, , u,- - sin2u7A;, 

$ + 3 ^^)$ + [(2cos2/x7fc- l)A; + 2( ^-^Y + 2^ [k cos 2 fijk ^-^) 

2/^7 2yU7 2/i7 

_-d sm^ = 47rG'(^5v9 - pK^(^)5(p) • (48) 
op fXJ 

Subsequently, we will consider the diffeomorphism constraint equation. The perturbed 
diffeomophism constraint with holonomy corrections is 



D[N'] = ^ / d'x5N^[pd,{6t5K',)-p{dk5K'J-{—^)6';{d,5Ei)+8nGndM- (49) 

The diffeomorphism constraint equation can be obtained by varying the diffeomorphism 
constraint with respect to the shift perturbation: 

8^G^^ = pdc{5tSK',) - p{d,5K',) - C^^!^)6l{d,5Ei) + SttCtt^^v^ = . (50) 

Using the expressions 5Ef ( |T3ll . ( H2ll and equation ( l30ll . the above equation reduces 

to 

^el^ + (^^^?^)$] = AnG^d^d^ . (51) 
2//7 

Finally, we will derive the Hamiltonian constraint equation. We note that after the 
variation with respect to the background lapse iV, the constraint equation will be second- 
order and can be neglected. So one can obtain the Hamiltonian constraint equation by only 
varying the perturbed lapse 5N 

^H^[N] 1 r , sin2/i7fc ^ ™ sin Ji-jk ^ 1 2 ^ ^.i.^n 

= . (52) 
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Substituting the expressions SEf ( IT3|l . SK^ ( H2i) and equation (130!) into the above equation 
yields the Hamilton constraint equation 
2^ o/Sin2/i7A; • ^ sin2/i7A;,2 ^,sinZ27A;,2 9 ,sinu7A:,2,^ . • r / xr i 

(53) 

In addition, using Eqs.( l43|) and ( l44l) . with the help of the background equations (129|) . 
( 130|1 and f l3T|l . the perturbed Klein-Gordon equation can be expressed as 

+ 2( )Sip - V25(^ + pV^^Sip + 2pV^^ - 4(/p$ = . (54) 

2/i7 

Now, we replace ^^^2^^^ Hubble parameter H in the perturbation equations ( 153|) . ( HHj) 
and ( 15T|) such that these equations can be reexpressed as 

V2$ _ 3H$ - - aC^^^Y + p|,(^^^^)2]$ = 47rG((^5'v9 + py^(<^)M , (55) 

$ + 3H<i> + [2e - ^ + 2tf - pA(!l!L^)2i$ = ATiG{05ip - pVd'^)5ip) , (56) 

op fi^y 

da{^ + H<l>) = 4nG^5^^a ■ (57) 

As we have emphasized in the introduction, among the three classical perturbations 
equations ([3]), (jlj) and ([5]) only two are independent. However, when the gauge-fixing has 
been done before deriving equations of motion in the presence of quantum corrections, it may 
not produce all terms correctly such that this consistency can not be maintained. In order to 
preserve the consistency for quantum corrected equations fl55l) . fl56l) and fl57j) . the additional 
correction terms must be required. The simplest way is only to modify the equations fl56p 
as follow by introducing some additional correction terms, 

i^+{3m+hM-i-p-^C^^^r]}^+[2M'+^^ = A7rG{^5ip-pV^{^)5v) 

M Op /i7 op /i7 

(58) 

The proof of consistency of these equations (l55l) . (!58ll and (157|) has been given in the 
appendix. Obviously, in the classical limit, /x — )■ 0, the equations fl55|) . fl58|) and f l57|l reduce 
to the classical cosmological perturbations (E]), (jl]) and respectively. In addition, we 
must also point out that since the gauge-fixing has been done before deriving equations of 
motion, the introduce of the additional correction terms is not unique. In order to obtain the 
more complete and unambiguous quantum corrected cosmological perturbations equation, 
we must consider the gauge invariant variables and derive these perturbations equations in a 
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gauge invariant manner, which is under progress. Combing these equations, one can obtain 
the following second order differential equation for $ 

(p M op fi'j 

+ [8tf + m-2h- aC^^Y - 2pU^^r + 2p-^e]$ = . (59) 

fl'J op flj ^ 

In addition, using the relation between the extrinsic curvature k and the conformal Hubble 
parameter HI ( I35p . one can obtain 

sin;U7^ 2 _ 1 ± a/1 - 4(/i7)2e2 



2(//7)2 



(60) 



If we denote Shi = - (™^)2 = - , which results from the holonomy 

corrections in the presence of the metric perturbation, and Sh2 = EI — A; — p ^ ( '^"^//'^ ) ^ , which 
be introduced by the requirement of the consistency, then the above second order differential 
equation can be further rewritten as 

$ _ + [6H + 2p'^^^ + + 2[2e2 + e + p^^H - 5,2 + 25^]$ = . (61) 

Up to now, we have completed the derivation of the cosmological perturbation equations 
in the effective loop quantum cosmology with holonomy corrections. 



V. DISCUSSION 



The effects of quantum gravity on structure formation, generally called trans-Planckian 



issues, have been investigated intensively (for example, we can refer to 260. In loop quantum 



cosmology, the analogous issues have also been investigated in Ref. 



ff 16 



16l-l20l|. However, in 



Ref. [19|, they assume that after a super- inflation phase, the universe underwent a normal 
inflation stage. Then they find that the loop quantum effects can hardly lead to any imprint 
in the primordial power spectrum. Although in Ref.jo] the scale invariant spectrum was 
obtained and the holonomy effects also leave their imprint on the power spectrum, only 



the holonomy effects from a fixed background were taken into account. In Ref. [20|], the 
cosmological perturbation equations with holonomy corrections were derived in longitudinal 
gauge. But they consider only the cases of large scale metric perturbations. In this paper, 
along the Hamiltonian approach we have derived the cosmological perturbation equation 
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for scalar modes in longitudinal gauge in the presence of holonomy corrections. In the 
presence of metric perturbation, we find that holonomy effects influence both background 
and perturbations, which contribute the non-trivial terms Shi and Sh2- Therefore, the 
holonomy effects will affect the power spectrum such that it is possible that the quantum 
gravity effects will leave their imprint on the cosmic microwave background observed today. 
In the future work, we will investigate analytically and numerically the characters of power 
spectrum in the presence of holonomy corrections, which might open a window to test the 
loop quantum gravity effects. 

In addition, when ignoring the additional corrections term Sh2, which be introduced by 
the requirement of the consistency, the second order differential equation become 

$ _ + 2(H - + 2[H - + + 2Shi]'l> = . (62) 

if LP 

here we have used the Klein-Gordon equation 

In this case, we can furthermore introduce the Mukhanov- Sasaki variable v = Then 
the cosmological perturbation equation ( !62i) reduces to 

v-V^v+ [(?)• - (2)2 + H - + 4Shi]v = , (63) 

ip if 

In momentum space, the cosmological perturbation equation 0631) can be written as 

v-[k^- AShi - mljf]v = , (64) 

where k denotes the momentum and fnljrj = (|)' — (|)2+EI— H^. Therefore, the cosmological 
perturbation equation ( l64|) can be effectively viewed as imposing such a modified dispersion 
relation at quantum gravity phenomenological level. Obviously, in such a modified dispersion 
relation, both the energy and momentum are bounded. Here, we point out that, in Ref. 
{2?! . Y. Ling et. al have also proposed a bounded modified dispersion relation, motivated 
by the isotropic homogenous effective loop quantum cosmology with holonomy corrections. 
Although both are bounded, they are also very different, implying we can not simply use 
the background corrections instead of perturbation corrections. In the future work, we will 
furthermore discuss the implications of such two modified dispersion relations. 

Our present paper is the first step towards studying the holonomy corrected cosmological 
perturbation equations in the presence of metric perturbation. Since constraints are modi- 
fied, the form of gauge invariant variables should change as well. Therefore, it is necessary 
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to study the perturbations with different gauges or in a gauge invariant manner in this 
formahsm, which is under progress. 
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Appendix A: The proof of consistency of these equations ()55|) . (|58|) and (|57|) 

In this appendix, we will give a proof of the consistency of these equations (ISSjl . ( 155]) and 
(!57|) . Without loss of generality, we will only derive the Eq. (158|) from the Eqs. (155|) and 
fl57|) . From the perturbation equation fl57|) . by taking the spatial derivation, we can obtain 

— (V^<^)+M\/^(^-A7vG^VHip = 0. (Al) 
ar] 

In addition, using the corrected Raychaudhuri equation ( 133|) and the perturbation equa- 
tion fl55]l . the term can be expressed as 

= 3H$ + [6H - 87rG((^2 + pV{^))] $ + 4TTG{^6ip + pV^{^)6ip) . (A2) 

Therefore, we can obtain the following expressions: 

A(v2$) = 3e<i> + [3e + 6e2-87rG'((^2^py ((/?))]$ 

drj 

+[i2ee - 8'kG{2^^ + pv{^) + p^v^i^))]^ 

+A7rG[ip6ip + ip6ip + pV^{ip)6ip + pV^{ip)6ip + pipV^^{ip)6ip] , (A3) 

ev2$ = 3e^<i> + [6tf - 87rGH((^2 + pV{^))]^ + 4ttGM[^^ + pV^{^)6ip] . (A4) 
In addition, we can also expressed the term AnG(^V^6if as following with the help of f l54p 
4:7TG^V^6ip = 4:TtG^[6^ + 2mip + pV^^{ip)6ip + 2pV^{ip)^ - 4(^<i>] . (A5) 
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Collecting all the above expressions (lA3p . (lA4p and f lA5|) . we can obtain the perturbation 
equation fl58|) by straightly calculating. Similarly, we can also derive equation f l55|) or ( 157|) 
from the remaining two equations. Therefore, among these equations above only two of 
them are independent. 
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